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Abstract 

A chiral random matrix model with locality is constructed and ex- 
amined. The Nielsen-Ninomiya no-go theorem is circumvented by the 
use of a generally applicable modified Dirac operator which respects the 
Ginsparg- Wilson relation. We observe the expected universal behaviour 
of the eigenvalue density in the microscopic limit. 

Introduction 

The success of chiral random matrix models, xRMM, in describing certain 
aspects of non-pertubative QCD indicates that the correlators in the Dirac spec- 
trum are universal and suggests that the nature of the spontaneous breaking of 
chiral symmetry does not depend on the detailed dynamical properties of the 
gauge field. (See Jackson and Verbaarschot Q and Wettig, Schafer, and H.A. 
Weidenmiiller For a comprehensive review, see Verbaarschot Q.) Here, 
we want to extend these considerations by constructing a local xRMM through 
the explicit inclusion of the Euclidean derivative terms of the free Dirac opera- 
tor. There are several reasons for the construction of such models. In practice, 
the range of validity of random matrix theory in describing real lattice gauge 
simulations is determined by the ability of the gauge field to mix free quark 
eigenstates. (In condensed matter physics, this range is known as the "Thouless 
energy".) The completely democratic treatment of all basis states in usual chi- 
ral random matrix models makes it impossible to address this question. For the 
same reason, local quantities, e.g., propagators, cannot be considered in usual 
xRMM. This limitation can be overcome by the construction of a local xRMM. 
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In order to represent the Euclidean derivative by a finite size matrix, it 
is necessary to discretize a finite (Euclidean) space-time volume. This leads 
inevitably to the problem of unintended fermion doubling. If one insists that 
chiral symmetry maintain its canonical form, {_D,75} = 0, on the lattice, this 
doubling cannot be avoided. In order to obtain a sensible local xRMM, we must 
resolve the fermion doubling in a manner which preserves a natural extension 
of chiral symmetry for finite volumes and lattice spacings. The impossibility of 
constructing a lattice Dirac operator, D, with locality and exact chiral symmetry 
but free of fermion doubling is a consequence of the Nielsen-Ninomiya no-go 
theorem Some years ago, Ginsparg and Wilson suggested the possibility 
of circumventing this no-go theorem by modifying the chiral symmetry condition 
to the form {D,75} = aDj^D. Following the recognition by Hasenfratz Q| that 
fixed point actions of QCD satisfy the Ginsparg- Wilson relation, attention has 
again been drawn to lattice theories which satisfy the Ginsparg- Wilson relation, 
see Hasenfratz, Laliena, and Niedermayer 0, Hasenfratz Neuberger [l0| , 
Liischer and Narayanan Liischer has shown that the Ginsparg- 
Wilson relation ensures an exact "lattice chiral symmetry" of the fermion action 
for any finite lattice size. Because of this symmetry and the fact that the 
explicit breaking of chiral symmetry enters through a simple condition on D, 
the Ginsparg- Wilson approach suggests itself as a suitable way of avoiding the 
no-go theorem in a local xRMM. Adopting this approach, it is possible to find 
a resolution of the fermion doubling problem and embed it in a local xRMM. 

The object of the present letter is to present a general procedure for the 
construction of a Dirac operator which satisfies the Ginsparg- Wilson relation 
and solves the fermion doubling problem. This procedure is applicable to lat- 
tice QCD as well as to the local xRMM which we construct. Further, we 
shall consider the spectral properties of this local xRMM and demonstrate the 
universality of the microscopic spectral density, i.e., the spectral density near 
eigenvalue zero on the scale of the lowest Dirac eigenvalue. 

The Ginsparg- Wilson relation and its implementation 

The Ginsparg- Wilson (GW) relation is a constraint on the Dirac operator D 



where a is the lattice spacing. The Dirac operator appears as the heart of the 
fermionic sandwich, ijjDijj, in the fermionic part of the QCD action 



Since the anti-commutator of D and 75 is not zero, chiral symmetry is explicitly 
broken. Instead, as shown by Liischer a new "lattice chiral symmetry" is 
present. The infinitesimal variation of the fermion fields associated with this 
new symmetry is 



Dj5 +15D^ aD-f^D , 



(1) 




(2) 



X 



ijj V + £75(1 - 7^ 



aD)il) and V^''/' + £'0(l 0^)75 • 
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Extension of this singlet flavour lattice chiral symmetry to a non-singlet lat- 
tice chiral symmetry is straightforward. Both variations are symmetries of the 
fermionic action provided that the Dirac operator satisfies the GW relation. 
Thus, the GW relation ensures the presence of a continuous symmetry which 
can be regarded as the remnant of chiral symmetry on finite lattices. The 
Nielsen- Ninomiy a no-go theorem states that lattice Dirac operator cannot si- 
multaneously satisfy a series of physically reasonable demands. Among these is 
that {-D,75} = 0. By explicitly modifying the form of chiral symmetry through 
the GW relation to the form of eqn. the no-go theorem can be circumvented. 
It is necessary to adopt a strategy for dealing with the fermion doubling problem. 
Here, we choose to work with the Wilson lattice realization of the continuum 
operator, H — ■y^ll ~ 7/^(9^ + but this choice is not mandatory. Let us 

first see how to implement the GW relation. 

Consider the Dirac operator, D, in an arbitrary matrix representation. Since 
D is to be multiplied by the row, tp, from the left and by the column, ip, from 
the right, D must be a square matrix of dimension, say, 2N x 2N. For some 
2N X 2N matrix e(-ff), we can write D as 

D = -ii-TAH)), (4) 
a 

wher^ 75 — diag(l, 1, .., 1, —1, —1, .., —1). Since 75 is invertible, there is a one- 



to-one correspondence between D and e(iJ). As noted by Narayanan |12 , the 
requirement that D satisfies the GW relation is equivalent to 

m))' - 1 • (5) 

Since the dynamics of D enters through e{H), one must build into e{H) some 
regularization of the nai've discretization which effectively eliminates fermion 
doubling. Having chosen a suitable form for H, we define the matrix e{H) as 

€{H) = Uh diag(sign(/ii), sign(/i2), .., sign(/i2Af)) Ujj (6) 

where Uh is the unitary matrix that diagonalizes H, 

uIhUh ^ di£ig{hi,h2,..,h2N) . (7) 

While any choice of H which provides a suitable regularization of the nai've 
discretization can be adopted, it is essential for the construction of e{H) that 
H be hermitian (or antihermitian) and that det(_ff) 7^ 0. 

Before turning to the construction of a local chiral random matrix model, 
it is useful to consider some general properties of the spectrum of the resulting 
Dirac operator and of the order parameter for the finite lattice chiral symmetry. 
Since e{H) is hermitian and satisfies eqn. (^, the combination ^^e{H) is unitary. 

^Hasenfratz |; | has shown that the index theorem is reproduced for finite lattice spacing 
and size provided that the GW relation holds. Thus, the difference in the number of 1 and 
— 1 entries in 75 depends on the topological index, u = — n_, which we shall take as zero. 
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This implies that D is normal and that the spectrum of D hes in the complex 
plane on the circle, (1 — e'^)/a with 9 G [—tt,7t]. Furthermore, D satisfies the 
hermiticity relation, ^^Dj^ — . As noted in this property in conjunction 
with the GW relation implies a complex conjugation symmetry in the spectrum 
oiD: 

D^x = X^Px (8) 
Djsi'x = A*75V^A ifA*^A (9) 
75Va = ±^A ifA*=A. (10) 

Using this complex conjugation symmetry, a Banks-Casher-like relation |0 
appears in the limits volume oo then m — > 0, where m is a regulator mass. 

(V^V) = {D-\x,x)) (11) 

1 r^" pis) , ^ 

' ' ds (12) 



m^O Vol^oo a^O Vol J j ^ (1 " 6°") ja^m 

1 Z""" 2mp(s) , 
= hm — / —J^ds , 13) 
m^o Vol Jq s^ + m^ ^ ^ 

where the last equality relies on the identity p(s) — p(— s), which is a conse- 
quence of the complex conjugation symmetry. The average in eqn. ( |lll) is over 
gauge field configurations. Finally, we obtain the Banks-Casher relation for the 
order parameter in the limit of zero lattice spacing 

m - ^ . (14) 

It is also straightforward to evaluate the integral in eqn. ( p^ ) in the limit 
m — > for any finite a. One obtains 

Hasenfratz ||^ has suggested the "subtracted" chiral condensate as a possible 
order parameter for lattice chiral symmetry, 

(^V').ub-^(Tr(i?-i-ial)) , (16) 

where the trace extends over colour, fiavour, and Dirac space indices. By virtue 
of the GW relation, only zero modes of D contribute to ('0V')sub. The result of 
eqn. (|l^) is that a Banks-Casher relation remains valid at all finite lattice spac- 
ings provided one uses the subtracted chiral condensate as the order parameter 
for the lattice chiral symmetry. Of course, in the limits of infinite volume and 
a — > 0, the lattice chiral symmetry becomes genuine chiral symmetry, and the 
subtracted chiral condensate evolves into the true chiral condensate. 
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A local chiral random matrix model 



We shall construct a local xR-MM using the standard Wilson prescription 
for dealing with the fermion doubling problem. We choose as a basis the states 
obtained by chiral projection from the eigenstates of the nai've lattice realization 
of the free Euclidian continuum operator 7^9^. In this basis, H is defined as 



H = 75X, where X = l 



A' iA 
iA A' 



iW 
iW^ 



(17) 



In general, the matrix W is uniquely determined by the gauge field configuration. 
The local xRMM Dirac operator is defined through eqns. (Q) and (||) by replacing 
W by a suitable random matrix. The N x N matrices A and A' are real and 
diagonal with 

1 ^ 

A = -diag(...,[^sin2(app^)+sin2(a^T)]i...) (18) 

^ i=l 
1 ^ 

A' = -diag(. ..,^(l-cos(app^)) + l-cos(a7rT),...) , (19) 



where —n/2 < j < n/2 — 1 with n = N^^^ . The spatial momenta 
pf'^ — 27rj/(an) are determined by imposing periodic boundary conditions in 
the spatial directions. Due to the anti-periodic boundary conditions in the 
time direction, the temporal momenta are given by the Matsubara frequencies 
{2k + l)7rT', where k is an integer. As is customary in chiral random matrix 
models, we have retained only the lowest Matsubara frequency, ttT. 

It is conventional in xRMM to replace both the known form of the free 
Dirac operator as well as the gluon field contributions in the random matrix 
W . The new feature of the present model is that the free Dirac operator is 
explicitly retained and W is assumed to describe only gluonic contributions. 
We choose the N x N matrix to be a random complex matrix, and its 
entries are drawn at random on a Gaussian distribution. This choice of W is 
motivated by the conjecture that the correlators of QCD with three colours 
are given by the chiral Gaussian unitary ensemble. (The Gaussian orthogonal 
and symplectic ensembles are expected to be relevant for the description of QCD 
with a smaller number of colours.) Note that the form of the M^-dependent 
part of X in eqn. ( p7| ) provides the most general antihermitian matrix which 
preserves the relation, 75X75 = X'^ , and thus ensures the hemiticity of H. The 
only "dynamical" information in W is its variance. We shall consider this point 
below. 

To see that this implementation of the GW relation yields a physically sen- 
sible non-doubled spectrum for the Dirac operator, it is useful to look at the 
free spectrum. Turning off the gluon field, i.e., W = 0, the eigenvalues of the 
free Dirac operator can be found analytically. Diagonalizing D, one finds that 
the eigenvalues, A, of D are in one-to-one correspondence with the eigenvalues. 
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S,j± — 1 ~ o,Ajj — aiAjj, of the free Wilson operator X. The result is simply 
A = (1 — £,/\£,\)/a- Evidently, only the phase information of ^ is carried through 
to A. This information, however, is all that is needed in order to resolve the 
fermion doubling problem: For small momenta, the Wilson term makes a negli- 
gible contribution to the real part of ^. The corresponding A eigenvalues appear 
at small angles, 9, on the circle (1 — e*^)/a. As the momentum increases, the 
Wilson term becomes increasingly important with the effect that the doubled 
states are pushed towards 9 = ±tt. Since the doubled states are thus well sep- 
arated from the small 9 region of physical interest, this represents a sensible 
resolution of the doubling problem expected to remain valid when W is reintro- 
duced. When W is not equal to zero, X is no longer normal and such analytic 
evaluation is not possible to our knowledge. (Of course, H = 75X is hermitian 
for both interacting and non- interacting fermions.) 

In general the analytical relation between the eigenvalues of D and the eigen- 
values of X, A = (1 — ^/\£,\)/cL, holds if X is normal and satisfies 75X75 — X^. 
Thus, this relation also applies when there is no deterministic part in X, i.e., 
A = A' = 0. The eigenvalues of D can then be expressed in terms of the 
eigenvalues, iwj, of the matrix 



iW 

iW'^ 



(20) 



which has the usual form of chiral random matrices. The eigenvalues iujj are 
simply mapped onto the circle (1 — e'^^)/a according to cos9j — 1/ + ojj. 
This is sufficient to ensure that both the microscopic spectral density and the 
various spectral correlators of D are identical to those of the original xRMM. 

Applications and numerical results 

The simplest application of the random matrix model constructed in the 
preceding section is the study of the chiral condensate as a function of temper- 
ature. As indicated by eqn. (^, this requires us to consider the microscopic 
spectral density of D. We wish to consider the case where chiral symmetry is 
spontaneously broken at T = 0. In ordinary random matrix models, the scale of 
the problem is set solely by the variance of the random matrix W. Chiral sym- 
metry is always spontaneously broken at T = 0, and the critical temperature for 
its restoration scales strictly with the variance of W. The situation is somewhat 
more complicated in the present case where D contains both deterministic and 
random elements. Here, the existence of a chiral condensate even at zero tem- 
perature is determined by a competition between low-energy free eigenvalues 
and the variance of the random matrix W. The spectral properties of matrices 
which are the sum of deterministic and random parts have been considered in 
refs. |, |l3 and (18). This work shows that the resulting spectral correlators 
and the microscopic spectral density will be given exactly by the strict random 
matrix results provided only that the random part of the matrix is sufficiently 
strong. We expect to find a similar result in the present case. 
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The problem considered explicitly in |^ was a matrix of the form 

A + M/t ) 

where A is a fixed diagonal matrix and the elements of W are drawn at random 
on the Gaussian weight 

P(W) ~ expi-NT.^TrWW''} . (22) 

The parameter E introduced here determines the variance of the distribution. 
When the diagonal elements of A are all non-zero, the chiral condensate, S, can 
be expressed as 

S = E^x , (23) 
where x is the only real and positive solution of 

iVE2 = y ^„ ^ , . (24) 

In the absence of a solution to eqn. (^J), the chiral condensate is zero. Noting 
that the right side of eqn. (|2j) is a monotonically decreasing function of x, 
it follows that a necessary condition for a non-zero chiral condensate in the 
presence of the deterministic matrix A is 

^ E • (25) 

Of course, neither the matrices X, i7, or _D of the present model have the 
form of eqn. Nevertheless, we would like to use eqns. ( p^ and ( psj ) to 

provide a rough estimate of the variance of W required to ensure a reasonable 
critical temperature for chiral symmetry restoration. These equations are nat- 
urally dominated by the lowest eigenvalues of A. Turning to the matrix X 
defined by eqns. we see that these lowest states correspond to the 

"undoubled" fermion states. Restricting our attention only to these states, it is 
reasonable to neglect the Wilson term. A'. This leaves us with a matrix having 
the form of eqn. (pi]). 

Having first determined the size, iV, of the matrix and the lattice spacing, 
a, we pick the variance of W in order to fix the estimated critical temperature 
at some value, T^. Thus, 

N/9, ^ 

Numerical simulations indicate that this procedure is sound in practice and that 
the resulting critical temperature is close to this estimated value. 

As an initial application of this local xRMM, we have performed a numerical 
study of the microscopic level density. As indicated above, we anticipate that 
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Figure 1: Universal behaviour of the eigenvalue density for the Dirac operator 
of the local xRMM. The histogram shows the microscopic level density at zero 
temperature while the smooth dotted curve is the analytic Bessel function ex- 
pression suitably scaled. This simulation was performed for a 24 x 24 lattice 
with a lattice spacing of 1/4 fm and a critical temperature of approximately 150 
MeV. 



the microscopic spectral density for our model will be identical to the usual 
xRMM result, which can be expressed in terms of the first two Bessel functions 
[|9| and [|0| as: 

p{X) = 2NE{NEX)[J^{2NXE) + JI{2NXE)] . (27) 

The analogous result for the Gaussian symplectic ensemble has been shown 
to be in excellent agreement with the results of QCD lattice simulations at 
low temperature. Eqn. (|2^) applies whenever chiral symmetry is spontaneously 
broken, and the magnitude of the chiral condensate is directly related to the 
asymptotic value of the microscopic spectral density. As anticipated, numerical 
results obtained for a variety of temperatures below Tc are in agreement with 
eqn. (|2^). Results for the p at T = are shown in the figure. The density is 
plotted as a function of 6 with A = (1 — ex-p[i0])/a. The disagreement between 
analytic and numerical results seen for larger values of is a consequence of 
both finite size effects and, more importantly, the limited ability of the random 
matrix, W, to mix free quark states of very different momenta. 
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Conclusions and discussion 

We have suggested a simple and well-defined procedure for modifying the 
usual Wilson Dirac operator in order to satisfy the Ginsparg- Wilson relation. 
The result is a lattice Dirac operator which both solves the fermion doubling 
problem and retains a lattice chiral symmetry for all lattice spacings and sizes. 
As noted, this lattice chiral symmetry becomes genuine chiral symmetry in the 
continuum limit. While we have used this procedure to construct a local chiral 
random matrix model, we emphasise that it is more generally applicable. (For 
example, the procedure outlined here would be suitable for QCD lattice simula- 
tions.) As is common in xRMM, we have introduced temperature dependence 
through the lowest Matsubara frequency. We have studied the microscopic level 
density of the resulting Dirac operator numerically and confirmed the antici- 
pated agreement with the universal behaviour found empirically in QCD lattice 
data and analytically in other chiral random matrix models. This allows us to 
extract the continuum limit (i.e., iV — > oo) of the lattice chiral condensate from 
finite-size simulations. The existence of this lattice chiral symmetry for every 
a leads us to anticipate that the limit a — s- 0, in which lattice chiral symmetry 
becomes genuine chiral symmetry, will be relatively smooth. 

This local xRMM represents an appealing intermediate step between pure 
random matrix models and real QCD lattice simulations. There appear to be 
a variety of interesting applications of local xRMM. The presence of the free 
Dirac operator (in a form free of fermion doubling problems) makes it sensible 
to consider a number of local properties which are genuinely not accessible to 
ordinary random matrix theory. As noted, it is possible to explore the interplay 
between random and deterministic elements in this model and to study the 
range of applicability of random matrix techniques to lattice simulations. It 
is now possible to search for quasi-universal behaviour in correlators and to 
address questions regarding localisation. While such studies will initially be 
numerical, they are far simpler than full QCD lattice simulations and may well 
have interesting insights to offer. 
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